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Abstract 

We consider the effect of interactions on tlie line shape of the two-photon 
Is — 2s transition in a (doubly) spin-polarized atomic hydrogen gas in terms of 
the interatomic interaction potentials. We show that the frequency-weighted 
sum rule for the intensity of the line is not given simply in terms of the 
pseudopotentials that describe the interactions between low-energy atoms. 
The origin of the departures from the simple pseudopotential result for the 
frequency-weighted sum rule is traced to what we refer to as incoherent contri- 
butions to the spectral weight. These arise from more complicated final states 
of the many-body system than the ones usually considered. In particular, we 
show how the relevant response function may be treated in a manner similar 
to the density-density response function for Fermi liquids, and express it as a 
coherent part coming from single particle- hole pairs, and an incoherent part 
coming from other excitations. We argue that in experiments only the coher- 
ent part of the response of the system is observed, and its contribution to the 
frequency-weighted sum rule is shown to be given correctly by the pseudopo- 
tential approximation. Finally we calculate the width of the coherent part of 
the line due to collisional damping. 

PACS numbers: 03.75.Fi, 32.70.Jz, 32.80.Cy, 67.65.-hz 
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I. INTRODUCTION 



After two decades of concentrated effort Fried et al. recently succeeded in realizing 
Bose- Einstein condensation in spin-polarized atomic hydrogen in a magnetic trap U|. In 
this experiment a key role is played by collisional frequency shifts, since the density of the 
atomic hydrogen cloud is monitored by observing the frequency shift of the Doppler-free 
peak in the two-photon Is — 2s absorption spectrum 0. Previously, collisional frequency 
shifts have also been observed in hydrogen masers and in atomic fountains [Q, where 
they lead to a serious limitation on the stability of these devices. A thorough understanding 
of such shifts is therefore central to the interpretation of various experimental results with 
atomic quantum gases. 

In the theory of line shifts that is currently standard one considers only the normal 
state of the gas and uses a Boltzmann equation to determine the effect of collisions on the 
absorption profile. The line shift is then found to be proportional to the difference of the 
Is — 2s scattering length, ais_2s, and the Is — Is one, ais-is- Recently two papers have 
appeared on the theory of the line shifts, one employing the random phase approximation 
0, and the other using sum- rule arguments |]^. One striking prediction of these calculations 
is that the line shift in a dilute, fully Bose-Einstein condensed gas should be one half that for 
an uncondensed gas of the same density. These papers have in common the assumption that 
the interaction between atoms may be assumed to be of the usual contact pseudopotential 
form, and that the interactions may be taken into account in a mean-field approach. In this 
paper we investigate the problem allowing for a more general interaction. We demonstrate 
that the frequency-weighted sum rule is given in terms of the hare interaction potential, not 
the pseudopotential. By using microscopic many-body theory we trace the origin of the 
discrepancy between the true frequency-weighted sum rule and the one calculated using the 
pseudopotential to incoherent contributions to the atomic propagators, which arise when 
a Is atom is excited close to another such atom. The latter processes, while relatively 
infrequent in a low-density gas, give contributions to the spectral weight at frequencies very 
different from those for excitation of an atom far away from any other atom. However, as 
we demonstrate in this paper, the shift of the coherent contribution to the response, which 
corresponds physically to excitation of an atom when it is relatively far away from other 
atoms, is given by the pseudopotential result. 

We have organized the paper as follows. In Sec. we first derive an exact sum rule 
for the frequency- weighted spectral weight, and will show that this is not given correctly 
by the pseudopotential result. In Sec. fT| we then study the problem from a microscopic 
point of view, and indicate how the absorption spectrum can be separated into coherent 
and incoherent parts. We also argue that the coherent part of the response is of greatest 
interest experimentally. In Sec. |V| we determine the collisional broadening of the coherent 
absorption peak and we end in Sec. |V| with our conclusions. 

II. SUM RULE APPROACH 

Let us begin by considering a system of hydrogen atoms in the Is ground state. The effect 
of applying the laser radiation is to excite some hydrogen atoms to the metastable 2s state, 
which has a radiative lifetime l/r2s of the order of one second. Experimentally, the hydrogen 
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clouds investigated are inhomogeneous, but since the length scale for density variations is 
large compared with the microscopic lengths in the problem, it is an excellent approximation 
to take the effects of inhomogeneity into account in the local density approximation, and 
consequently in our calculations we consider a spatially uniform system. If the radiation 
field is spatially uniform, its interaction with the hydrogen gas may be represented by a 
perturbing Hamiltonian 

ifi = ^ I rfx [e--VL(x)^i.(x) + e^-VL(x)^2s(x)] , (1) 

where u is the angular frequency of the pair of photons, the operators V^^lx) and '?/'a(x) 
create and destroy atoms in the state and is the effective Rabi frequency determined 
by the strength of the laser field and atomic matrix elements. 

The unperturbed part Hq of the Hamiltonian is given by the sum of the intrinsic atomic 
energies of isolated atoms at rest, the kinetic energy associated with the translation of atoms, 
and terms that take into account interactions between atoms. To an excellent approximation 
the interaction energy is given in terms of local two-body potentials dependent only on the 
distance r between atoms, and we denote the potential for two atoms in the Is state by 
Vis-isir) and that for one atom in the Is state and the other in the 2s state by Vis-2sij). 
Since we consider the case of weak excitation, we shall not need to specify the interaction 
between two excited-state atoms. In detail we thus have 

Ho = J d:>i V'L(x) f — ^ + ei. J ^is(x) + j dx iplsi^) f — ^ + 62, J V2.(x) 
+ ^ J d:s.J dx.' tpl,{:>i)tpl^{yi')Vis-is{:>^ - x')V^is(x')^/'i^(x) 

+ JdxJ rfx' ^L(x)V^L(x')Vi,_2.(x - x')^2s(x')^i.(x) , (2) 

where m is the mass of an atom and denotes the energy of the atomic state a. Note that 
for clarity we have in the first instance neglected the effect of the finite lifetime of the excited 



atom. In Sec. |IT4 however, we show how it can be easily incorporated into the theory. 

The net rate of transitions may now be calculated from Fermi's Golden Rule, and is 
given by 

27r 

= YT.\HHi\n)\^6{hij + En - E^){pn - Pm) . (3) 



m,n 



Here Pn is the initial probability for occurrence of the many-body state |n), which is an 
eigenstate of the Hamiltonian Hq and therefore obeys HqIu) = En\n). We note that for the 
situations of interest in the Bose- Einstein condensation experiments, initial states containing 
2s atoms play essentially no role, since the probability of 2s atoms being present is very small 
because the energy difference between a 2s atom and a ground-state one is much larger than 
the thermal energy fc^T. From Eq. (|^) we thus find that the rate of absorption of energy is 

271 

htoliuj) = — Y.iEm - En)\{m\Hi\n)\''5{huj -Em + , (4) 

m,n 

and the average frequency of the line is given by 
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_ _ / dujul{u) _ Y.m,n{Era - En)\{m\Hi\n)\^Pn , . 

Jdujl{uj) flEm,n\{MHi\n)\'^Pn 

To evaluate the average frequency we, following the procedure adopted by Oktel et al. 
0, again make use of HQ\n) = En\n), and consider the thermal average of Hi[Hq, Hi], or 
equivalently the double commutator [Hi, [Ho,Hi]]. In contrast to Ref. 0, we however do 
not assume that the interaction may be represented by a pseudopotential. The average 
frequency is then given by 

_ {Hi[H„Hi]) 

^ = n{H!) • 

The physical content of this equation is that the average frequency shift is given by the 
difference in energies of the expectation value of the energy in the initial state and that in 
the state created by operating with Hi on the initial state. Evaluating the expectation value 
of the commutator expression above directly, we find for the frequency shift relative to its 
value for an isolated atom the result 

T) f 

ALU = -jdr [Vu-2sir) - ri,_i,(r)](72(r) , (7) 
where n is the density of the gas and 



\ (4sir)4smumisir)) (8) 



is the pair distribution function for ground state atoms in the initial state of the system. In 
arriving at this expression we have again neglected the possibility of 2s atoms being present 
in the initial state. This result is simple to understand, since the operator Hi converts a 
single ground-state atom in the initial state into an excited state one with an amplitude that 
does not depend on position. The average energy difference between the initial state and 
the one created by the laser is therefore the energy required to convert a Is atom into a 2s 
one. Since the masses of the atoms in the two states are the same, there is no contribution 
from the kinetic energy, and the sole contribution, apart from the energy difference for an 
isolated atom, comes from interactions. This situation should be contrasted with that of 
an isotopic impurity, like a ^He atom in liquid "^He, which is just the opposite, in that the 
masses are different, while the interaction potentials are the same. The sum rule derived 
here is analogous to sum rules for spin response of condensed matter systems, and for spin, 
isospin, and spin-isospin response of nuclei. In these cases the basic origin of the shifts is 
terms in the interaction that are not invariant under rotations in spin, and/or isospin space, 
or, in the present problem, rotations in the pseudospin space corresponding to conversion of 
a Is atom into a 2s one. 

The long-wavelength assumption is appropriate for the two-photon transition when the 
two photons that are absorbed have equal and opposite momenta. When the total momen- 
tum q of the absorbed photons is non-zero, the perturbing Hamiltonian depends on space 
and we need to generalize the sum rule to spatially varying interactions. This is straight- 
forward and we find that the average frequency shift is given by adding the recoil energy 
h^(^ /2m to the spatially homogeneous result for q = 0. 
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Let us now compare our result in Eq. (0) with that of earher work. If the interaction 
potentials are weak, the correlation function will vary little over the ranges of the potentials, 
and we may replace the pair distribution function by its value for zero separation. We then 
obtain 

(A^)b = -^?2(0) j dr [V,s-2s{r) - Vu-is{r)] . (9) 

This is equivalent to the result of Levitov et al. [0], since for weak potentials the Born ap- 
proximation may be apphed, and thus the scattering lengths a are related to the interaction 
potentials by A-Kfi^a/m = J dr V{r). 

The interaction potentials for hydrogen atoms are not weak, and the Born approxima- 
tion is not valid. Therefore it is important to explore how the pair distribution function 
behaves at short distances. On length scales larger than the range of the atomic interac- 
tions, correlations should be well described in terms of mean fields. However, for strong 
potentials it is not permissible to assume that the correlation function for small separations 
varies slowly on distances of the order of the range of the potential. Rather one expects that 
the many-body wave function for small particle separations will behave as that for a pair 
of atoms interacting via the Is — Is interaction, since the effects of other atoms will then 
be negligible at low densities. Hence, provided the energies of elementary excitations of the 
system are small compared with the typical energy scale over which the two-atom relative 
wave function changes significantly, it will be a good approximation to assume that the pair 
distribution function scales as the square of the scattering wave function at zero energy, i.e., 
l^reK'")!^- We therefore write 

^2(r)^|vl/,ei(r)rr7r(0), (10) 

where the mean-field correlation function g^^{0) is the pair correlation function on length 
scales that are greater than the range of the interaction but small compared with other 
lengths in the problem, such as the thermal de Broglie wavelength, the particle separation, 
and, when a condensate is present, the coherence length. We have chosen the normalization 
of the wave function such that at distances large compared with the range of the Is — Is 
potential it behaves as "^reiir) ~ l — au-is/r. Note that the above procedure is equivalent to 
assuming a wave function of the Jastrow form to describe the correlations at short distances, 
the Jastrow factor being taken to be of the form of the relative wave function of two atoms 
at zero energy [Q. 

The final result for the shift is thus 

= ^92''{0) J dr [Vu-2s{r) - Vu-is{r)]\%ei{r)\' . (11) 

This expression cannot be simply rewritten in terms of scattering lengths. A simple example 
that demonstrates this is a Is — Is interaction with a hard core at a radius Tc- The expression 
for the frequency shift does not depend on the Is — 2s potential at distances less than r^, 
since the relative wave function for two atoms in the |ls) state vanishes there. However, the 
Is — 2s scattering length is sensitive to the behavior of the Is — 2s potential at distances 
less than Tc, and therefore this is incompatible with the frequency shift being expressible 
solely in terms of scattering lengths. We expect the Jastrow form of the wave function to 
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be accurate irrespective of whether or not the gas is Bose condensed, and consequently in 
a completely Bose-condensed gas the shifts are predicted to be a factor of 2 smaller than 
in a gas of the same density with no condensate, reflecting the usual 2! reduction factor for 
two-body processes §]. 

The result of this calculation is that the frequency-weighted sum rule is quite different 
from what one predicts if one uses the pseudopotential. To understand the origin of these 
differences it is convenient to explore the problem from a microscopic viewpoint. 



III. MICROSCOPIC APPROACH 

To understand the sum- rule result, it is helpful to think about the nature of the final 
states that can be created from the initial state by the operator of interest, which in this case 
converts a Is atom into a 2s one. Relative to the initial state, the simplest excited states 
have an extra 2s quasiparticle and an extra Is quasihole, and will be referred to as single 
quasiparticle-quasihole pair excitations. In the random phase approximation these are the 
only states take into account. The physics of the process may be understood by regarding 
the degree of freedom associated with converting a Is atom into a 2s one as a pseudospin. 
If the commutator of the pseudospin-raising operator with the unperturbed Hamiltonian is 
zero, there is a unique frequency for all transitions. Because the interaction between a Is 
atom and a 2s atom differs from that between two Is atoms, however, the Hamiltonian is 
not invariant under rotations in pseudospin space, and its commutator with the pseudospin- 
raising operator is not zero. Consequently there can be transitions to states with a range 
of energies. It is perhaps helpful to consider a spin system in an applied magnetic field. If 
the interaction between the particles commutes with the spin-raising operator, the raising 
operator will couple only to states whose energy differs from that of the original state 
by h times the Larmor frequency. However, if the interaction is not invariant under spin 
rotations, other excited states with different energies can be created. In Fermi-liquid theory 
the first sort of transitions correspond to the creation of a single quasiparticle-quasihole 
pair, while the more complicated excitations correspond to creation of many pairs. For the 
problem under study here, an n-pair excitation has one extra 2s quasiparticle, n extra Is 
quasiparticles, and n + 1 extra Is quasiholes. The difference between the results for the 
frequency-weighted sum rule calculated with the pseudopotential and the actual potential 
is due to the multipair excitations. For Fermi liquids an analysis of the density response in 
terms of single-pair and multipair states may be found in Ref. |]T0[|. A formulation of the 



problem for more general sorts of response was presented in terms of microscopic theory by 
Leggett , and the results were discussed in terms of Fermi-liquid theory in Ref. . 



Let us begin by expressing the result for the transition rate in terms of the response 
function for the operator 

O = i I dx (e-'-VL(x)^i.(x) + e^-VL(x)^2s(x)) , (12) 

where V is the volume of the system. The response function is defined in the usual way as 
the temporal Fourier transform of the retarded commutator, and is given by 

XM = di^j^ dt e^-*([^L(x,t)^2.(x,t), 4.(0, 0)^1,(0,0)]) 
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where p„ is again the probabihty of the state n being occupied and we neglected the occu- 
pancy of the final state compared with that of the initial one. The transition rate in Eq. 
is therefore given by 

m = -l"^ [n(a^)] , (14) 

where 

liioo) = ^^\{uj) (15) 

is the polarizability of the gas. This is the desired result, because it explicitly shows that 
the transition rate is related to the polarizability of the gas, which is easily accessible with 
equilibrium many-body techniques. Indeed, in that language hli{uj) is equal to the (retarded) 
self energy for the "effective photon" causing the Is — 2s transition and the imaginary part 
therefore determines its finite lifetime, which physically is due to absorption by Is atoms in 
the gas. We are thus left with the task of calculating the polarizability, which theoretically 
implies that we have to evaluate the diagram in Fig. |I[ We begin by considering two simple 
calculations, the Hartree-Fock approximation and the random phase approximation, before 
discussing the more general formulation. 



A. The Hartree-Fock and Random Phase Approximations 

In this section we consider a number of examples where only coherent contributions to 
the response are taken into account. These calculations lead to results identical with those 
of Oktel and Levitov P]. To familiarize ourselves with the present formulation, let us first 
consider the ideal Bose gas. Then Eq. (|13]) becomes 

X{uj) = ~J d^j^dt e^-*(V^L(x,t)^i,(0,0))(^2.(x,t)V'L(0,0)) . (16) 
Moreover, the single-particle propagator is given by 

(V'L(x, t)V>i,(0, ^)) = Y.Y. iVke-^'^-^+^(^''+^--^)*/'* , (17) 

^ k 

where ek = h^\<^ /2m is the kinetic energy of a Is atom, yU is the chemical potential for 
Is atoms, and A^k = 1/ [e^^'^'^~^^'> — 1) is the Bose distribution function with (3 = l/fc^T. 
Similarly we have, including now the finite atomic lifetime of the 2s atom, 

(V^2.(x, t)V^L(0, 0)) = :^ E e*'^-"-'(^''+^2^-'^^^=/2-/.)tA . (18) 

^ k 



Substituting the last two results, we find for the polarizability in Eq. ([T5|) the expression 
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Therefore, we conclude that the absorption hne of the gas has a profile given by 

Kuj) - iV^^^)' ^ (20) 

which is just the number of atoms = nV times the atomic line profile and exactly centered 
at the atomic resonance in this case. Note that diagrammatically we have now calculated 
the lowest order contribution to the polarizability in Fig. |1|, in which the exact Is and 2s 
propagators are replaced by the ideal gas ones and there are no vertex corrections. 

At the next level of approximation we dress the Is and 2s propagators by including the 
effect of atom-atom ladder diagrams as shown in Fig ^ This corresponds to a Hartree-Fock 
approximation, in which the effective interaction is taken to be the T matrix for two-body 
scattering. For a gas with no condensate, the effect of dressing the propagators in the above 
calculation is to replace by ea + ^S^, where to lowest order in the T matrix 

^S,. = !!L^l£zi£^ , (21) 
m 

and 

hE,^ = l![^i£z^ . (22) 
m 

The factor-of-two difference between the numerical factors in Eqs. (^Tj) and ( p2D reflects 
the fact that both the Hartree and Fock terms contribute to the energy of a Is atom, but 
only the Hartree term contributes for a pair of unlike atoms. Because these interaction 
corrections to the atomic energies are purely real, the absorption line is of the same shape 
as in Eq. (^) , but is now centered at a frequency shifted from the single-atom resonance by 
an amount 

(Aa;)HF = {ais-2s - 2ai^_is) . (23) 

m 

This is the "naive" Hartree-Fock result for the coUisional frequency shift due to the mean- 
field interaction that a Is and a 2s atom experience from the surrounding gas of Is atoms. 
Most important for our purposes is that if we repeat the above calculation for a fully Bose 
condensed gas of Is atoms, we find that now 

(Ac^)hf = (ais-2s - ais-is) , (24) 

m 

which shows that only the contribution from the Is — Is mean-field interaction is reduced 
by a factor of 2, in agreement with the fact that Is and 2s atoms are distinguishable. 

We have called the above Hartree-Fock approximation naive, because it is well known 
that for an approximation to satisfy the conservation laws, it is necessary to include vertex 
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corrections in addition to the self energy corrections discussed above. For this problem this 
amounts to including the effects of the mean field self-consistently. Formally, these two kinds 
of corrections are related by the condition that the vertex correction must be the functional 



derivative of the self energy corrections with respect to the applied field ||T3[ . In our case this 
implies that we also have to calculate the "maximally crossed" diagrams shown in Fig. |^. 
These correspond to the chains of particle-hole bubble diagrams calculated in the random 
phase approximation and in Fermi liquid theory. This is easily achieved since it corresponds 
to summing the geometric series 



huj - (Au;)hf + eis - (^2s + ifi^2s/'2 

1 47rais-2s^ n 1 

hu - {Auj)}iF + (^is - <^2s + iKT 2s/ 2 m hu - {Auj)}if + (^u - <^2s + iKT 2s/ 2 

In the end we thus find that 



+ 



^^""^ ^\\hu;- (A^)rpa - {e2s - e,,))^ + {hT2s/2y ' ^^^^ 



with 



(AcuJrpa = (At^jHF H = [ais-2s - ais-is) ■ (26) 

m m 

Moreover, in the fully Bose-Einstein condensed case the vertex corrections are absent and 
we recover the Hartree-Fock result, 

(Au;)rpa = (Aa;)HF = {ais-2s - ais-u) ■ (27) 

m 

We therefore conclude that for a fully condensed Bose gas the collisional frequency shift is 
indeed reduced by an overall factor of 2, in agreement with our sum-rule result in Sec. |T| 
and the work of Ref . . 



B. Coherent and Incoherent Contributions to Response Functions 

To explore the physics in greater detail, it is convenient to adopt an approach exploited 
in the context of Fermi-liquid theory [|1T|. We first expresses the single-particle propagator 



in a many-body system as the sum of a coherent part coming from an intermediate state 
with a single quasiparticle excitation, and an incoherent part coming from more complicated 
excitations. Mathematically this implies that 

G(2)(p, e) = G2t(p, 6) + gS(p, e) , (28) 

where p is the momentum, and e is the energy. The coherent part, which corresponds to 
the quasiparticle, is given by 

G'lM = ^^y (29) 
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where Z(p) is the renormahzation factor or quasiparticle residue, and e(p) is the quasi- 
particle energy. Both of these quantities depend on the atomic species considered. The 
incoherent contribution to the propagator corresponds to transient effects due to the dress- 
ing of a free atom to make it into a quasiparticle. Likewise the two-particle propagator 
G*^^^(p, e; p', e') for a pseudospin fluctuation may be expressed in terms of a coherent part, 
corresponding to a single quasiparticle-quasihole pair, plus an incoherent part coming from 
multipair excitations, i.e., 

G^'\p,e;p',e') = Ji't^, Ji^'^l,, + G^{p,e;^,^) . (30) 
(e-ei.(p))(e'-e2s(p')) 

We next analyse the diagrams for the response function x by separating the single particle 
propagators into their coherent and incoherent contributions, as was done by Leggett in the 
context of Fermi systems. We then divide these diagrams into two classes. The first class 
contains those diagrams which are reducible with respect to the coherent contributions of 
two (one Is and one 2s) single-particle propagators. We refer to these as the coherent 
contribution. The second class contains all diagrams which are not reducible in this sense 
and we call this the incoherent contribution. We remark that since the operator O does not 
change the total particle number, the two coherent particle lines must have their arrows in 
opposite directions, and therefore correspond to a quasiparticle-quasihole pair. 

Expressed in a formal matrix notation the response function may be written as 

X = Tr[G'(^)(l-r(4)G(^))-i] , (31) 

where T^^^ is the two-particle vertex function that is irreducible with respect to two particle 
lines with oppositely directed arrows, i.e., it is irreducible in the particle-hole channel. Sep- 
arating out the terms that contain only incoherent contributions to G^^^ from the others, we 
find 



X = Xinc + Xcoh , (32) 

where 

Xinc = tv[gS£(i - r^^^GfS)-^] (33) 

and 

xcoh = iv[(i - r(^)Gf2)-GSi(i - rfi gSI)-^(i - r^'^cl^r] ■ (34) 

The factor (1 — F'^^'^gI^^)"^ corresponds to a vertex renormahzation and the quantity 

rSl = r(^Hi-GS£r(^^)-^ (35) 

is a renormahzation of the interactions between the coherent parts of a particle-hole ex- 
citations due to intermediate states with incoherent particle-hole pairs. For the present 
problem, an important feature of this result is the existence of the incoherent contribution 
to X, since this is what is responsible for the difference between the sum rule evaluated with 
the pseudopotential and the true sum rule. 
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We turn now to the coherent contribution to the response function. The coherent part of 
G^^^ has the same form as for two particles with energies modified by the medium, apart from 
the renormahzation factors Z. However, if one multiphes the matrix element for coupling 
of the two photons to the excitations by a factor (1 — T^'^'^ G^^^)'^ {Z2sZisy^'^ and uses for 
the effective interaction between a quasiparticle and a quasihole the quantity ZisZ2sT^^\ the 
coherent contribution to the response has precisely the same form as in the random phase 
approximation calculation above. This modified interaction plays a role analogous to that 
of the quasiparticle-quasiparticle interaction introduced in Fermi-liquid theory by Landau. 

Let us now analyse the consequences of the above for a low density gas. In that case 
the renormahzation factors tend to unity, and the two-particle vertex reduces to the T 
matrix. The quasiparticle energies reduce to the Hartree-Fock ones, and the mean-field 
interaction is also just the T matrix. Thus the coherent contribution to the response has 
precisely the form predicted by the mean-field theory calculation in Sec. [Ill A| . Observe 



that in calculating the average frequency associated with the coherent part of the response, 
the renormahzation factor for the effective two-photon matrix element cancels out. What 
implications does our calculation have for experiment? In addition to a sharp peak in the 
absorption due to the excitation of a single quasiparticle-quasihole pair, the calculation 
predicts a broad background due to creation of more complicated final states. However, 
because the background is expected to be a rather smoothly-varying function of frequency 
that extends over a large frequency range, it is difficult to detect. Consequently, the part of 
the absorption spectrum that is investigated experimentally is only that due to the coherent 
contribution to the response function. 



IV. COLLISIONAL BROADENING 

We now consider how collisions broaden the coherent part of the line. One effect is 
that the self energy of the atoms acquires an imaginary part, and the coherent parts of the 
propagators become 

^^<P-)- e-e(p)S-(p.e(p)) - ^^'^ 

where S"(p, e) is the imaginary part of the self energy. Another effect is that there are vertex 
corrections analogous to those responsible for the contributions to the line shift beyond what 
is predicted by the Hartree-Fock approximation. The total width is most easily calculated 
by observing that the propagation of a 2s atom and a Is hole is determined by the difference 
between the Is — 2s interaction and the Is — Is interaction. For definiteness, let us consider 
a gas with no condensate. In the absence of the Is — 2s interaction, the only contribution to 
the width comes from the imaginary part of the self energy of the Is atom, which is given 
in the dilute limit by |T^ 

^KsiP, ep) = -2n ( — ^^^^^ ) H ^(^p + V " ^p+p" - V-p") 
\ / p'p" 

X [A^p,(l + Arp,„p„)(l + ATp+p,,) - (1 + A^pOArp,_p„Arp+p»] , (37) 
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the factor of two being the result of the Bose enhancement of the cross section, which is due 
to the exchange process. Since the imaginary part of the self energy is momentum dependent 
the absorption line of the gas is in principle not exactly Lorentzian. Nevertheless, the typical 
width of the line is determined by the average — (2/iV) J2p ^p^'LiPy ^p) ^^^1 thus equals 

~ nrn^ ^ ^\^p + V - ep+p" - ep'-P"J 

pp'p" 

X iVp[iVp,(l + iVp,_p„)(l + iVp+p.) - (1 + iVpOiVp,_p.iVp+p»] . (38) 
When the Is — 2s interaction is included, the result is simply 

Ai 2. - 2^ d[ep + Cp/ - Cp+p" - ep/_p"j 

pp'p" 

X iVp[iVp,(l + iVp._p.)(l + iVp+p") - (1 + Np')Np'^p"Np+p"] . (39) 
In the classical limit, this reduces to 

AT2s = 87m{ais-ls - ais-2.)^(frel) , (40) 

where (frei) = 4(A;T/7rm)^/^ is the average relative velocity between two Is atoms in the 
atomic hydrogen gas. The total width of the line is the sum of the natural width and the 
collisional contribution, and it is therefore equal to r2s+Ar2s. The width is of order {ais-2s — 
ais-is)/'^T times the shift, where At = h/{2nmkBTy^'^ is the thermal de Broglie wavelength. 
When a condensate is present, the above calculation can be easily generalized. In a first 
approximation we only need to take into account explicitly the macroscopic occupation of 
the zero-momentum state by substituting A^p ^c^^p.o + ^p, with ric the condensate 
density. At the next level of approximation we also need to incorporate the Bogoliubov 
coherence factors. 



V. CONCLUSIONS 

In this paper we have considered the effect of interactions on the two-photon absorption 
line profile in spin-polarized atomic hydrogen by means of a frequency-weighted sum rule and 
by means of microscopic many-body theory. We have shown that the line profile consists 
of a narrow coherent peak on top of a broad incoherent background. For typical atomic 
potentials this background in principle has sufficient spectral weight that the pseudopotential 
approximation does not give an accurate estimate of the total contribution to the frequency- 
weighted sum rule. However, the frequency of the narrow peak, which is the feature most 
easily seen experimentaly, may be expressed in terms of the low-energy pseudopotentials. 

We have also shown that the collisional frequency shift of the absorption line is reduced 
by a factor of two if the gas is fully Bose condensed. We have pointed out that for this 
factor-of-two reduction it is crucial to take many-body correlation effects into account that 
go beyond the Hartree-Fock approximation commonly used for these dilute atomic gases. 
At this point it is worth mentioning that the Bose-Einstein condensation experiments by 
Fried et al. apparently do not to see this effect. Their results seem to be consistent with 
the Hartree-Fock theory, which, due to the fact that ais-u <^ |ctis-2s| for atomic hydrogen, 
basically predicts no reduction at all At present we have no explanation for the cause 
of this discrepancy. 
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FIGURES 




FIG. 1. The polarization diagram that determines the two-photon absorption hne shape. The 
thick Unes denote the exact Is and 2s propagators, and the smah and large black areas denotes 
the bare and exact vertex functions, respectively. 




I f 

FIG. 2. The Is and 2s propagators in the T-matrix approximation. The thin lines represent 
the bare propagators and the dashed lines the interactions. 




FIG. 3. The vertex corrections in the T-matrix approximation. 
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